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Abstract 

We present a formalism able to generalise to a relativistically co- 
variant scheme the standard nuclear shell model. We show that, using 
some generalised nuclear Green's functions and their Lehmann repre- 
sentation we can define the relativistic equivalent of the non relativistic 
single particle wave function (not loosing, however, the physical con- 
tribution of other degrees of freedom, like mesons and antinucleons). 
It is shown that the mass operator associated to the nuclear Green's 
function can be approximated with the equivalent of a shell-model 
potential and that the corresponding "single particle wave functions" 
can be easily derived in a specified frame of reference and then boosted 
to any other system, thus fully restoring the Lorentz covariance. 



1 Introduction 

The difficulties one meets in building a theory for relativistic bound system 
with finite number of particles are well known. Up to now, in spite of many 
efforts in this field (see |2| for a comprehensive review, but also [3] to get 
an example of the present approach to the relativistic shell model), to rec- 
oncile relativity, translational invariance and shell model seems to be a very 
hard task. Moreover, even the connection between exclusive and inclusive 
processes is non-trivial for two order of problems. 
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By one side in fact in a relativistic framework it is impossible to fully 
disentangle the nucleonic and nuclear dynamics |3] (in a few word the nu- 
cleon form factors do not factorize) even in the simple scheme of the Plane 
Wave Impulse Approximation (PWIA), because the separation between lon- 
gitudinal and transverse motion is a frame-dependent concept and the Fermi 
motion of the nuclei prevents a full separation of them in a nuclear context. 
Moreover even the concept of Coulomb sum rule as it is usually interpreted 
looses its meaning and can be regained only at the prize of introducing a suit- 
able renormalization factor (that, fortunately, turns out to be largely model 
independent) [5j. 

On the other side, when going beyond PWIA multiple counting of dia- 
grams occurs jHj, with the consequence that the integral of the differential 
cross sections for (e, e'p) or (e, e'n) reactions no longer coincide with the 
inclusive cross section (this because of the existence of channels where, for 
instance, another nucleon is emitted but not revealed). 

Thus with the usual many-body expansion it is very difficult to connect 
forward scattering amplitude with the total cross section. This suggested us 
to extend the idea of Green's functions (of any kind) by allowing situations 
where the kinematics of the initial and final nucleus can be different. This 
will be suitable to directly study the elastic processes in a fully Poincare 
invariant way, but natural extensions could also be obtained (and we plan do 
pursue this line in the future) by choosing different initial and final states. 

For the moment we limit ourselves to the problem of two interacting 
particles, namely a nucleus ± a nucleon or (if case) an elementary particle 
(nucleon) ± a quark. This job enables us to account for recoil effects in 
high energy nuclear reactions and in quark physics (excitations of nucleon, 
meson) . 

We first begin in section |2] with a short review of what it happens in the 
non-relativistic frame, in order to provide a layout of the matter we would 
like to generalise, and also in order to make easier the understanding of the 
origin of some problem we are concerned with, i.e., if they arise from the 
many-body theory or from the relativity. 

Next, in Section |3] we consider the one particle (or hole) problem in pres- 
ence of a nucleus (to be more specific we consider nuclei with A nucleons ± 
one nucleon). We use a formalism similar to the one of the Green's function 
formalism. Due to nuclear recoil the equations for hole and particle become 
different unlike the infinite systems or models where recoil is neglected. 

In Section E] we introduce a model with particle (hole)-nucleus interaction, 
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conceptually similar to the shell model. The equation of motion in such 
approach looks like the equation for a particle embedded in a mean field, but 
the equation is relativistically covariant (of course with suitable choice of the 
interaction). The crucial point here is that we need to introduce the "shell" 
ground state of the system. 

Section provides a possible perturbative scheme to go beyond the mean 
field level, and sec. El presents a simple case where the "relativistic shell 
model" can be easily solved, thus coming in contact with the real world, not 
pursuing abstract and inapplicable theoretical formalisms. 

2 The non-relativistic single particle motion 
in nuclei 

The purpose of this section is twofold: by one side to remind the reader 
the general scheme of some old models for the propagation of a particle (or 
hole) inside a nucleus, that pertain to the history of nuclear physics, but 
that (partially) could be scarcely useful today in practical calculations; on 
the other side we wish to remind those topics that can provide a guideline 
for our generalisation to a relativistic finite nucleus and to "non-diagonal" 
Green's functions (this topic well be clarified in the next section). 

If we limit ourselves to the propagation of a particle (or hole) inside a 
nucleus the most natural framework to begin with is certainly the Feschbach's 
approach to the optical potential [7J Let us repeat once more its main 
topics (or, better, the ones we need in the following). The general idea is 
that, if a Hamiltonian is defined in a Hilbert space 7i, then we can project 
Schrodinger equation into a subspace 7i' C TC the price to be payed being an 
energy dependence in the effective potential. As everybody knows, if V is a 
projection operator 

V : TC — >TC' 

and Q = I — V then the Schrodinger equation in the space TC' reads 

>= EV\^ > 
(1) 

The conceptual points we want to remark are the following: 

• the optical potential can be defined in this way and its most relevant 



H opt (E)\q>-- 



VHV + VHQ 
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structures can be derived, but eq. ((JJ can by no means be used to 
evaluate it. 



• Eq. (JTJ) is quite general: according to the definition of V it can be 
adapted to a variety of problems: we shall consider in the following the 
particle and hole propagation in a nucleus but we shall also remind its 
application to (e, e'p) reactions in impulse approximation. 

• Eq. specifically imposes causality at each time. We shall see that 
this "microscopic" causality is the first reason that inhibits a micro- 
scopical calculation of the optical potential. 

• The optical potential displays an imaginary part, but since (JTJ) is de- 
rived from a true Schrodinger equation in a bigger space, the eigenvalues 
are necessarily real. This applies of course to the discrete ones, i.e., to 
stable nuclear states. In order to conserve a Lehmann representation 
with real eigenvalues, the usual way out is that of discretizing the whole 
system by means of a box normalisation. 

The next requirement is the definition of V, i.e., of the physical problem 
we will concerned with. In the archetypal case, namely the elastic scattering 
of protons off nuclei V reads 

V = J d 3 rd 3 r'p(r,r')^(r)\$ A >< ®° A \ip{r') (2) 

where \<& A > is the ground state of a nucleons with A nucleons, ip and ifi 
are the non-relativistic destruction and creation operator of a nucleon in the 
point r (spin and isospin will be neglected throughout this paper) and p is 
fixed by the requirement V = V 2 , that implies symbolically 

P = J J — (3) 
1 — n 

where 

n(r, r')=< ^^(r^r'M > . (4) 

The space of the solutions of the eigenvalue equation ((TJ) is of course 
isomorphic to L 2 (R 3 ). We can define an orthonormal basis in L 2 (R 3 ) by 
defining 

|r) = J <Pr' {I -n}~^ {Y,Y')i)\r')\§° A > . (5) 
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We also assume that this basis is complete. The eigenvalue equation now 
reads 

J d 3 r'{T\H opt {E)\v')ip n {r') = E<p n (r) (6) 

and we are in position to connect the solutions of © with the true eigenstates 
of the A + 1 system: let > a solution of the complete Schrodinger 

equation in the space of A + 1 particles. Some simple algebra enable us to 
write the solutions of © in the form 

<p n (r) = J d 3 r'{I-n}^ (r,r') < $>(r')|<n + i > (7) 

the eigenvalue being of course E n . For future reference let us define the 
function 

M')=<*A + i\tf(r)\*A> (8) 
We have shown above that up to a rescaling ip n is solution of the eigenvalue 
equation for the "optical" Hamiltonian. It also follows that ip n is, by itself, 
eigenstate of the (unsymmetric) operator 



Lopt 



Very much in the same way we can handle other situations. In particular 
we shall be concerned with hole propagation inside a nucleus. Thus we define 
another projection operator, namely 

V = J d 3 rd 3 r'ip(r)\$° A > rT^r') < {r') (9) 

and all the above formalism follows up provided we do the substitutions 

I — n — > n 

and 

Vn(r)— ^„(r)=<$:U|V(r)l<> ■ ( 10 ) 
The next step, in the early times of the optical potential, was the con- 
nection with the single particle Green's function of the system, defined, as 
usual, as 

r( n < ^ A \T{jj(x) } ^(x')}\^ A > 

G ( x ' x ) = Z ^oi^n ^ • ( n ) 



It was shown by Bell and Squires jU] that the self-energy (or mass operator) 
can be interpreted as an optical potential (not coincident however with the 
one introduced by Feshbach and discussed above). As is well known G can 
be separated into a retarded and an advanced part having the following 
Lehmann representation: 



G = G + + G- (12) 



G {t,t,u) - > _ (13) 



n 



frU*(r A 



n 



where the functions ij) n and n are those discussed above. Note that the 
functions ip n and <p n are in some way connected with the shell model: in 
fact the index n runs over all the possible nuclear eigenstates, but grouping 
together some levels and constructing in this way the single-particle levels 
and discarding those with a too small strength one is lead back to the shell 
model. This however implies the breaking of the translational invariance, 
since the latter would strictly imply the functional dependence 

n ,^n ~ e ip ' r . 

How to recover the translational invariance and at the same time to leave suf- 
ficient room t2o introduce the analogous of the "shell-model wave functions" 
will be the task pursued in the following. 

The previous discussion enables us to write down (but by no means to 
solve or to approximate) the inverse of G ± . We can write indeed 

[G ± ] _1 (r, r» = u - T - M ± (r, r'; uo) (15) 

where M.^ could be called the mass operator for particles or holes, with the 
property 

G+CW = G-ft(r) = (16) 

for uj = E^ ±l . The discussion above shows that ultimately 

G ± = VT^H opt ^=L= (17) 
V 1 — n 
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provided the particle or hole projection operator is used in the rhs. In this 
way we have indirectly defined the mass operators Ai ; it must be reminded 
however that the whole Green's function obeys the Dyson's equation 



where the mass operator (or self-energy) can be derived from a perturbative 
expansion, but it is not the sum of the two mass operators defined sepa- 
rately for particles and holes, and they can in no way be derived from any 
perturbative scheme. 

Before ending this section we would also remind that the same formalism 
have been employed in studying (e, e'p) reactions in the frame of the Distorted 
Wave Impulse Approximation (DWIA). There it is convenient to define many 
projection operators, any of them pertaining to a residual nucleus left in an 
excited state plus an outgoing nucleon. This again is formally correct but by 
no means one can be able to explicitly write down the (almost) infinite set of 
different optical potentials. Thus one ultimately ends up with assuming the 
same optical potential for the outgoing nucleon independently of the state 
the residual nucleus is left in. As a non-trivial consequence the differences 
between longitudinal and transverse channels are lost (see JU] to recover 
it). Again here the most natural treatment of the problem goes through the 
introduction of a particle-hole Green's function having the form (we follow 
the standard notations) 



where j M is the electro-magnetic current. The differences (we could better 
say the incompatibility) between this approach and the DWIA has already 
been shown by the authors of this paper in ref. [Bj. 

3 The generalised one-body Green's function 

In a relativistic approach, with the aim of pursuing the analogy with the de- 
scription of the non-relativistic single-particle or single-hole motion discussed 
above, and moreover in order to avoid the disease of multiple counting of dia- 
grams as outlined in [HJ , we consider a bound system of fermionic and bosonic 



G = G + G MG 



(18) 




(19) 
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fields with finite baryonic number and in the ground state in its frame of ref- 
erence but assuming in general non-zero different total momenta for incoming 
and outgoing states. 

This approach will be applied here to nuclear physics, but the study of 
the quarks dynamics in a nucleon could also be an affordable task. Moreover 
in both cases the accounting of the recoil effects is allowed. 

Let us first of all introduce the incoming and outgoing nuclear bound 
state \p > and p' > for a nucleus of mass M and initial and final 3-momenta 
p and p'. We can of course introduce initial and final 4-momenta by putting 
Po = a/p 2 + M 2 and p' = a/p' 2 + M 2 . The normalisation reads 

<p>\p>=(2n) 3 2p 5 3 (p-p') . (20) 

We now define a generalised single particle Green's function as 

G pp ,(y,y') = —^=<p'\T{^y)My>)} \p> . (21) 

We have already observed in the previous section that in a 2-points Green's 
function the translational invariance will rigidly constrain the analytical form 
of the functions n and ip n thus forbidding its interpretation (within some 
approximation schemes) as single-particle wave functions. The above choice 
of writing a generalised single particle Green's function, with, actually, one 
more argument, relaxes the above constraint and will turn out to be the key 
issue in constructing the relativistic generalisation of the shell model without 
violating the Poincare invariance. 

The particular case we are considering deserves a comment about the 
realization of the linked cluster theorem. To understand it we could interpret 
G pp i(y, y') as a limiting case of a two-particle Green's function. Imagine that 
^(p, t) is the destruction operator of the nucleus in its ground state with 
total momentum p and |0 > the physical vacuum. Then G pp i(y,y') can be 
regarded (up to normalising factors) as the following limit: 

< 0|T{^),^0,* t (p,0,*(P^')} lo > 
lim hm : 

t^-oot'^+oo < j > 

and the denominator < 1 > is the tool that ensures the cancellation of 
the disconnected diagrams in the two (composite) particle Green's function. 
Of course this property is preserved through the limiting process and conse- 
quently G pp '(y,y') (where we neglect the denominator < 1 > throughout 
this paper) has always to be intended as constructed by linked diagrams only. 
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Now we want to represent the function G pp >(y,y') in Lehmann represen- 
tation. First however we need some kinematical considerations. 
If p is the 4-momentum operator, i.e., p = {po = H, p), 



Pa\P >= Pa\P > 



then we know that 



i>{y) 



if>e' 



-ip-y 



(22) 



(with if) = ip(y = 0)). With these definitions we can write the Fourier 
transform of G pp '(y,y'). Here however an ambiguity arises, since our G is 
ultimately, as quoted above, a two-particle Green's function, we can choose 
in Fourier transform two different kinematics, one tailored for the propagation 
of a particle and one for a hole. For the "hole" channel, whose kinematics is 
depicted in fig. [I] we define 



p-q 



p -q 



p-q 
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p ' A p A p' A 

Figure 1: The kinematics for the "hole" propagation 



G h (p } p';q,q') 



d 4 ydye^-^-^'-^'G pp ,{y,y') 
{2nf5\q-q')G^p'-q) 



(23) 



Let us remark that all the formalism here is covariant, and in order remind 
this property we indicate in Gh a dependence upon the 4-vectors p, p' and 
q. Strictly speaking the mass shell condition on p and p' would entail a 
dependence upon the 3-vector part only 

The Lehmann representation of Gh is easily written down as 

<hM; «) = ^ < p¥#^W + f^^Ul* > (24) 

^VPoPo H-q -ia q — H + ia 

where 

q = p + p - q , (25) 

that reflects the kinematics of fig. H There the intermediate lines denote the 
sum over all the eigenstates of the system with baryonic number A — 1 (first 
term) or A + 1 (second term in the rhs of fig. having total 4-momentum 
q. The former states are characterised by 

H\q x >= q\o\q\ > ■ (26) 

where A is an index running over all the A — 1 states, their mass being M\, 
with the relations 



<?A = (<7A0,q), g A0 (q) = ± V / M A 2 + q 2 (27) 

the normalisation being 

<gv|gA>=(27r) 3 ||5 AV 5(q-q'). (28) 

The states with A+l particles are denoted with the index v and is understood 
to run over all the A + l excited states. For them eqs. (J26j) to (J28|) also hold 
up to the replacement A — > v. 

Here, in order to be completely covariant we have allowed the intermediate 
states to contain negative energy solutions too. This case will be practically 
irrelevant in nuclear physics, but not at all negligible if we want to extend 
this formalism to QCD. 

Note also that, just to make a choice, we have assumed for a state with 
baryon number A a boson normalisation (A is assumed to be even). Thus 
consequently an A — 1 state must be normalised as a fermion. 

In order to make the Lehmann representation for Gh more explicit we 
introduce in eq. the complete set \q\ >< q\\ (for the A — 1 system) 
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in the first term of its rhs and ^ v \q v >< q u \ in the second one. We easily 
finddUE! 

G,(P,M = E ^ ( P', q) ^ (P ' q) + E ^ (P ' q,fe( P-q'' . (29) 
, 9Ao(q) - Qo - ^ Po+Po-%- 9i/o(q ) + 

where the and ■?/> are defined as 



r„(p.q)) = J-^-<p^| gw > . (31) 
z PoQuo 



The equal time commutation relations imply 

-4= < p'I {^(y)} Ip > e-^-y = 

< p'Ip > / ^=r^) = ( 27r ) 3j (p " p') • ( 32 ) 

Inserting now a complete set of intermediate states ^2\ u \q\,v >< Qx,u\ in 
the two terms of the anti-commutator in the lhs and using (}2*2*|) we get the a 
completeness equation in the form 

E pa(p, <M(p'> q) + E ^(p'> ^(P, q') = (2^) 3 5(p - P') • (33) 

A 1/ 

We observe that now the functions <p\ and ip\ play the same role of 4> n 
and ip n in the eqs. (|13|) and (|14j) of sec. El but now the formalism is Poincare 
invariant and further, even if we are only considering the nucleonic Green's 
function, all the information about the dynamics of the system is already 
embedded in ip x and ip\. 

In the same line as above, we can also introduce a "particle" kinematics: 
in analogy with we introduce 

G P (p,p';q,qi) = jd%d 4 y'G PP iy,y')e liq - p ' )y - tiqi - p)y ' (34) 
= (27ry6 4 (q- qi )G p (p,p';q) 
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q-p' 




whose graphical representation is given in fig. |2j 

By comparison with eq. (f2*3*|) one immediately establishes the link 

G P (p,p;q) = G h (p,p';q) . q' = p + p' - q (35) 

For sake of simplicity we consider a specific model of a fermionic field 
interacting with scalar bosonic field <j(y), neglecting the self-interactions ~ a 3 
and a 4 . In this simplified scheme (that nevertheless still contains all the 
difficulties relevant to the fermionic sector) the Hamiltonian of the system 
reads 

H = J d 3 y^ y (-ijV y + m + ga y ) ip y + H° a (36) 

H® being the free Hamiltonian of the a meson). Using the equation of motion 
for the field operator ip 

7 °[V',fll = 7[V',P] + (m + <7)V (37) 
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we can derive the evolution equation for G, namely 

(17 -dy-m) G pp , (y, y') = (38) 
(2vr) 3 5(p - p')5\y - y') - i—^== < p'\T {a{y)^{y)W)} \p > 

We can first of all prove that, on general grounds, that G pp >(y,y') can 
be inverted and consequently a mass operator can be defined by means of 
a perturbation expansion. The standard proof requires to introduce the 
generating functional for connected diagrams and then to perform a Legendre 
transformation on it, and since it can be found in the usual textbooks [IB"] , 
is not reported here. We only remark that the generalisation of the Green's 
function definition used in the present paper will only affect the boundary 
conditions of the path integral representation of the generating functional, 
but not the steps needed to define the mass operator. 

We rewrite eq. (}2"3*j) in the form of a Dyson's-like equation in Fourier 
transform as 

{7 • (p - q) - m} G h (p, p'; q) = (27r) 3 5(p - p') (39) 

[Mb.(p,Pi] q)G h (p 1 ,p; q) + M p (pi,p'] g') G pb>Pi; <?')] , 

the mass operators M.^ and M. p being defined according to the "hole" and 
"particle" channels through 

/ ^%M h ip,p 1 ;q)G h {p 1 ,p';q) = < p'\<p J® = q) , a<p\p> 

(40) 

/ ^MMM)G P (j>^q') = ^L=<p'\ya <P\P> 
J (2tt) j ZVPoPa qo-H + ia 

(41) 

In the above p and p', as well as p±, are restricted to the mass shell and 
a is defined as 

" = *(l/)U <r(y) = e*»*e*» . (42) 

Further, the index h m M.^ only remind the "hole" kinematics chosen in 
introducing the Fourier transform. In the configuration space the mass op- 
erator is univoquely determined by the Green's function and embodies both 
particle and hole propagation. 
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Concerning the structure of the mass operator, the general theory tells 
us that it is built by the sum of all 1PI (one-particle-irreducible) diagrams 
defined in sec. 03 as shown in fig. El and must have the analytical structure 



p-q 

pA 



< 




< 


> = 




> = 



p -1 

P'A 



p-q 

pA 



< 




< 




S lw!e 




> 




> 



p-q 
P'A 





Figure 3: The diagrammatic representation of the Dyson's mass operator 



Mh(p,p',q) = S h (p,p',q) + yt Thn f 1 



Fhnip, q)Th* n (p', q) 



/a 



(43) 



Here the first term is a smooth (regular) function of go while the second term 
carries poles in q^, living of course in the region q 2 > M|_ x . Eq. (|3*§|) also 
shows that a pole of M.^ corresponds to a of and vice versa. 

The "particle" mass operator has an analogous expansion but its poles 
lie in the region q 2 > M\ +l . 

The knowledge of the Green's function (or of the mass operator) gives us 
access to many observables, like, for instance, the number of baryons 



A = Ti 



d 4 q 



■7°Gh(p,p;g)e" 



(44) 
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and the ground state energy 



< p\H\p > 



Po = 



< p\p > 



+ 




(45) 



Now we can introduce the analogous of the "single particle wave func- 
tions" in our relativistic approach, generalising what was described by eqs. 
()7I8|) . This is done by taking eq. (ffilj) and integrating it over q in a small 
circle containing a pole q\o (within box normalisation if needed) and remem- 
bering that the poles of the Green's function and of the mass operator never 
coincide. In so doing we immediately find 



The above quantities <p are not, of course, wave functions, because they 
feel the presence in the system of antiparticles as well as of mesons, but can 
be looked at as eigenfunction of the system. The case of uniformly invariant 
system (free Fermi gas or maybe quark-gluon plasma) may enable us to 
make strongly simplifying assumptions. For finite systems we however can 
still exploit the idea of a mean field calculation. 

4 The relativistic shell model 

The last equations of the previous section contain the ground idea to build 
the relativistic analogue of the shell model. 
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We first consider the "hole" channel and rewrite eq. (|46j) in the form ^3] 



{l(p- qt(q)) -m}^i(p,q) = j j^^M h (p,p'; q)<p t (p', q) (48) 

with the subtle difference that now go is considered a free parameter. If 
follows that (|4~Kj) considered at a given g can be regarded as an eigenvalue 
equation, the eigenvalue being qeo{q) that is in general different from the go 
fixed and contained in A4 h . Here notations matter: in fact if£ depends upon 
the 4-vector g, chosen by the exterior. We have left the dependence upon p 
instead of p to remind the reader that (pi is a 4-spinor depending, further- 
more, upon Lorentz-covariant quantities like g 2 and p ■ g, being understood 
however that po is fixed by the mass shell condition. 

Having distinguished between g and the eigenvalue qeo(q), (}4~%j) will have a 
complete orthogonal set of eigenfunction, i.e., the ipe must obey the properties 

d?p 

(p* t (p,q)(pe(p,q) = S e > e c e (q) (49) 



(2tt) 3 

J2-^w(p,qM(p',q) = (2tt) 3 <5(p-p'), (so) 

cg(q) being a suitable normalisation factor. 

Eq. (J4*H|) at a fixed and suitably chosen g looks like a shell model equation 
having a (non-local) "shell model potential" Aih] the functions <pe(p, q) are 
not connected with any observable quantities, but they are expected, for a 
reasonable approximation of M.\y and in a convenient range of go (it means 
some average of the single particle levels of a shell model well) to approach 
the "single hole" wave functions <pe(p, q) previously introduced. This is of 
course likely below the Fermi level. 

For the particle channel we rewrite eq. (|47|) as 

f d 3 p' 

{l(Qi(q) -p)-m} xi(p, q) = / — .„ ; — =M p (p', p; q\)xe(p', q) , (51) 

J 2(2tt)VPoPq 

normalisation and completeness relations being fully analogous the "hole 
wave function" case. 

According to the general form (}2~9*|) we can now introduce the shell-model- 
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q'=p+p'—q 



like Green's function as 

ct{q) qeo{q) -q -iat 

(52) 

where q and de are suitable normalisation factors and 

N e = 6(£ F - £) , (53) 

with 

Y,N e = A. (54) 

i 

Of course the labels £'s are ordered increasingly with the corresponding en- 
ergy and £f denotes the highest occupied level (Fermi level). 

The equation the function Q^ m '(p,p'] q) fulfils is similar to (j3HJ), provided 
the 5 function in the r.h.s. is replaced by 

Yl ^° i —T\W(P' ?)^(p'» + / ,s Xt(p', q')xt(p, q') 

The above is of course expected to approximate a 5 as far as the "shell model" 
approximation holds valid. 

In the above formalism clearly the poles of £?hoie must also be poles of 
Q^ m ' (the converse is not true in general because some poles of are killed 
by the projection operator 1 — Ne). 

We assume that the equations 

Qa(QXt) = 9A,o(q) (55) 

Qtofavt) = ^o(q) (56) 

have one or more (maybe infinite) roots for a given "single particle" quantum 
number £. The ?a«o(q) are the eigenvalues of the equation for the A— 1 particle 

state 

H\q\e,L>e >= q\ e v e o\q\ e ,is e > ■ (57) 

The residua of and of Q^ m ' below the Fermi level coincide and from eq. 
and (J52J) we obtain 

Ne l\ dqeoiq)}^ 1 



^(p,q)^(p'>q) = j 1 - } Mp>q)n{p',q) 



(58) 
9o=<?Ao(q) 
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Thus the above suggests to attribute to the ip\ e the meaning of a generali- 
sation at the relativist level of a single particle wave function in a shell model, 
fully maintaining, nevertheless, Lorentz and Poincare invariance. We re- 
peat once more that this occurs because we are considering a "non-diagonal" 
single-particle Green's function where the recoil of the daughter nucleus is 
accounted for. 

Now we make a physical assumption that further narrows us to the shell 
model: we assume that a one particle level I is composed by the same sub- 
levels of the exact many-body problem in such a way that 



(2vr) ; 



= N t ' 

<?o=9A„o(q) 



This property is not peculiar of a relativistic system, since the same will 
happen in the non-relativistic case. 

Now we are ready to make the last step and introduce a phenomenological 
"shell model" potential V(p, p'\ q), with p and p' restricted to the nucleus mass 
shell and the 4-vector q chosen as 



q = (?o, q) , go = yq 2 + m\_ x , (60) 

where Ma-i is the mass of the daughter nucleus in its ground state (this 
choice maintains the analogy with the non-relativistic case: see, e.g., [T3] ) 

Thus the potential V is independent from q and we assume it to be sym- 
metric, i.e., V(p,p'; q) = V(p',p; q). From now on we must guess V(p,p'; q) 
on phenomenological grounds in such a way that 



f d 3 p' 

{lip -Qt(q)) -m}ipt(p,q) - J j^^V(p,p';q)ip e (p',q 

d 3 p' 



d 3 p' 

■ \ 

(2tt) 3 [ Mh ( p ' p> i q } ~ ViP'P'i^WiP'iV) 

(61) 



will be reasonably small. 
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Once a parameterisation for V has been given we can write down the 
eigenvalue equation for the "hole wave functions" 



{7 • (p - &(q)) - m} (ft{p, q) = / — — ^(p^'; q)^(p', q) (62) 




being 



e« = (fto(q), q) , p = (po(p), p) , p' = (pd(p')> p') • ( 63 ) 



The index £ (of course discrete) summarises now all the quantum numbers 
pertaining to a given "one-hole" state in the "shell model potential" V. 

As an aside, in analogy with the "hole" channel, we can introduce an 
equation for the "particle" channel as 



Up to now covariance has been preserved. The next step is to find a 
practical way to solve the "shell model equation" (|62|1. Since it is not so 
easy to give an explicit solution of it in a covariant form, we are forced, 
in the following, to choose a suitable reference frame where the equation is 
particularly simple. Once the solution has been found, however, we need a 
procedure to boost it to any reference frame. This will be our next task. 

Thus, coming to the "hole" channel, there are two natural choice for the 
reference frame. One is to assume it as the rest frame for the A — 1 system, 
i.e., q = 0: the eigenvalue equation there becomes 



h = (po(k) - q m , k) , k = (p (k), k) , k' = (p (k'), k') . (66) 



Another customary choice is to assume p = (rest frame for the A-nucleons 
system). Of course any reference frame can be reached by means of a boost. 
Thus let |0 > be the p = frame and let A be the boost from |0 > to the 
state \k > corresponding to q = 0: 



{7 ■ (Qe(<l) -p)-m} xe(p, q) = / -pr^V{p,p'; q)Xi{p', q) • (64) 





(65) 



where 



k >= A|0 > 
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Let us specify in more details the state of the daughter nucleus: we expand 
the previous index t as 

i = (A, J, M, tt) 

where J is the total angular momentum, M its third component and it the 
parity. A will then resume all the other intrinsic (not frame-dependent) quan- 
tum numbers. We can now write 

<Px,j,mAK 0) = J— = < 0; A, J, M, tt|^(0)A|0 > 
V 2 Po(k) 

= J_ < 0; A, J, M, tt|AA- V(0)A|0 > (67) 
V 2 Po(k) 

= J_ = S(-v) < -77 A k; A, J, M, tt|^(0)|0 > 
V 2 Po(k) 

where as usual S is defined through the relation 

S(-v)^(0) = A"V(0))A (68) 

and reads 

Sf v) J pM + M ( 7°7-k \ 

Of course 

Po( k ) 

denotes the velocity of the boost from the rest frame of the daughter nucleus 
and an extra factor accounting for the mass difference between the A and 
A — 1 systems is required, namely 

Vx = ■ (71) 

The above entails 

A _1 |0; A, J, M, tt >= | — rjxk; A, J,M,n> . (72) 
In the p = frame the 4-vector q transforms into 



Po(k) = yJiMffi + ^ k 2 (73) 



q = -- =^ = -ifck (74) 



20 



and the "eigenfunction" <p reads 



= 5(-v A )v9 Ai j iAf>7r (-k/r7 A ,0) 

with 

V A = """TTT • (76) 

(ZAJM7r;0W 

This definitions implies 

d 3 k 



/d k 
j^VxjmA ' - k )^A,j,A/,^(0, k) 



/d 3 k 
7^p^A,j,M,^(k, 0)v3A,J,M,^(k, 0) . (77) 

The formalism above has shown that we can solve the "shell model" 
equation (Jo^j) in the rest frame for the A — 1 daughter nucleus and then, using 
the above kinematics, transfer the solutions to the usual rest frame of the A 
nucleus, namely p = 0. Having established good transformation properties 
of the solutions, we now need to find them in the preferred reference system. 
Before doing explicit (model) calculations let us investigate a little what lies 
beyond the shell model. 



5 The perturbative expansion 



The shell model in nuclear physics is usually thought as the th order (mean 
field) of a perturbation expansion. Using the eigenfunctions derived from 
eqs. ()62|) and (JMJ) we can represent the "unperturbed" Green's function in 
the "hole" kinematic, in analogy with (J52|) . as 



Q(<?) *o(q) -qo-ia 
1 - N £ xe(p', <l')Xt(p, q') 



(78) 



q'=p+p'—q 
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Also, we can introduce a "shell-model particle Green's function" by ap- 
plying the relation to eq. (J7SJ), namely 

Gl(p,p';q) = G° h (p,p';p + p'-q) . (79) 

The first order in the expansion of the mass operator A4 in terms of 
the meson interaction coincides the Hartree-Fock approximation: the first 
contribution reads 



(A^h)Hartree = -tTrg 2 J jL^G h ° (p, p' , k)V (p - p') 

v— \ f d 3 k 



(80) 



where of course 



v o{k) = T2 — rr - (81) 

k — mr a + icx 

is the free a propagator, and is displayed in fig. HI while the second term 




G h o(p, P ';k) 

Figure 4: Diagrammatic description of the Hartree contribution to the mass 
operator 

represents, as one can easily convince himself, the Fock contribution, namely 
(M h ) Fock = iTrg 2 J -0^G h °(p,p',k)V o (k-q) . (82) 
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G h o(p,p';k) 



Figure 5: Diagrammatic description of the Fock contribution to the mass 
operator 



and is represented in fig. The other terms have complicated and in practice 
non manageable expressions that involve the detailed structure (i.e., the ex- 
cited states) of the target nucleus. We only show diagrammatically a second 
order contribution in fig. El 



6 A simple model 

The above theory looks rather formal. Thus let us show how it can be 
implemented in a practical case. In order to have manageable formulas we 
consider a "shell model potential" of the separable form 

V(k,k>) = -(27r) 3 ^y^,(k)y j t i (k')^(|k|)^(|k , |) (83) 
jij 

where as usual 

1 1 

= < l > J ' ~ s; 2' S '^' 2' J ' J > Yl i- s ^ Xs 

s " " 

are the generalised spherical harmonics and v j is some function to be chosen 
in such a way to reproduce the nuclear phenomenology. Actually we put 

" j(|k|) = ( iOTfr+igv) )' < 84) 
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where b an c are constant and a(A, J) will depend upon the atomic number 
A and, a priori, upon the total angular momentum J. 

Now we can solve eq. (jo^j) (we assume that, as established in sec. HI 
once the eigenf unctions ip n have been found in the frame of reference q = 
then the above described transformations can provide their expression in any 
other frame). 

The spinor (not strictly speaking a wave function) solution of (J65|) will 
be labelled by J and u and has the form 

^=MM («) 



Yji.j(k)G Jfl (k) 



where 
and 



, r . / M J 1 • ( 86 ) 



+1 for states with parity (— 1) J+7 
-1 for states with parity (—1' 

Using the well known relation (independent of the parity) 

<T-kY Jl3 (k) = -kY Jl , j (k) (87) 

we find 

, ^f[(ko-m)F J ^k) + kGj,Uk)}Y Jlj (k)\ 
( ' 7 ^ ~ V- [kFjAk) + (ko + m)GUk)] Y JVj {k)) m 

and now not only eq. (|fi5j) is separable thanks to the choice of the potential 
but the large and small components of the spinors are decoupled and one 

gets 

(h - m)Fj^(k) + kGj^k) = -vj{k)Xj[Fj^] (89a) 
-kFj^k)-(k + m)G.^(k) = -Vj(k)Xj[GjJ (89b) 

where the functional Xj is defined as 

oo 

Xj[f] = [ P 2 dpvj(p)f(p) . (90) 
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Figure 6: A typical example of higher order diagram in a perturbative cal- 
culation 
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Then inverting the above and recalling (J66|) . i.e., expliciting the expression 
of the eigenvalue k = po(k) — qm, we get the system 

t? iu\ {p (k) - q eo + m)Xj[F JyUJ } + kXj[Gj, w ] 

Fj^k) = -— ^ vj(k) , 91a 

(po(k) - q eo y -k 2 -m 2 

r kXj[F J)U) ] + (pp(k) - qep - m)Xj[Gj^} 

(Po(k) - qtoY - k 2 -m 2 

Inserting then these expressions into the definition of the functionals Xj[Fj ;UJ ] 
and Xj[Gj ;UJ } we get for them a homogeneous system, namely 

Xj[Fj^} = Bj{q lQ )Xj[Fj^-Cj(q m )Xj[Gj^, (92a) 
X.j\G, huj ] = Cj(q m )Xj[Fj, u ]-B+{q m )Xj[Gj, u ] (92b) 



where we have defined 

R ±/ \ 7,2,, [g<o-Po(k)±m]^(fc) 

5jM = /fc ^(po(k)- %0 ) 2 -fc 2 -m 2 ' (93) 



o 

oo 



CM.) - /^ ^-gU-^ ' (94) 



The eigenvalue equation generated from the system (|9*2*|) is then 

^jfeo) = [Cj(qm)} 2 - {Bj{q m ) - l)(B+(q m ) + 1) = . (95) 

Note that the functions and Cj are real only in the range qgo > Ma + m 
or qiQ < Ma — m, (the former case referring to a A-nucleus system plus 
an antinucleon and the latter to a A-nucleus plus a hole). Here of course 
Ma = Po(0) denotes the rest mass of the A-nucleus. 

Once the equation is solved in q^ we also get, up to a normalisation 
constant, the explicit expressions for the "wave functions" 

F,Jk) = - gffl - (^±^£M -*d zSM) „; (fc) , (96a) 

(po(k) - %o) -k 2 -m 2 

Gjjk) = - (to - ■»)(! - 5M) + *g&B> »» . (96b) 

(po(k) - % ) -k 2 -m 2 

Note that, from ([95)1. the energy levels are degenerate with respect to j and 
to the parity and in the notations we can rewrite qeo as qjo- 
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If we further introduce the notation 



e = M A - q J0 - m (97) 
then the functions B and C become 



B+(e) = - I k 2 dk (K ^ ■ ^ UA AP+ /, U v 2 j{k) (98) 



o 

oo 



(Ap + e)(Ap + e + 2m) - k 2 



Ap + e + 2m 2 ( 



B~(e) = - k 2 dk— -f- ■ , 7^vi(k) (99) 

' (Ap + e) (Ap + e + 2m) - k 2 jy ' y ' 



2, 



and 



C(e)) = / k 2 dk— — r^viik) (100) 

K " J (Ap + e)( y Ap + e + 2m)-k 2 JK J v ; 

o 

(101) 

where we put 

Ap = p {k)-M A (102) 

to better control the orders-of-magnitude: this last quantity is in fact ex- 
pected to be small (say, of the order of k 2 /2Ma) unless we look to extreme 
situations, and for bound states e is of the order of few MeV. 

To exemplify how the above works, we have chosen the parameters in 
PI) as 

1 r 12 
a(A,J) = — 0.7314 + 0.3274^3 - 0.0884s + 0.0089A - 0.005(2J - I) 
m . 

(103a) 

b = 0.09 (103b) 
c = 0.1 (103c) 

and we have evaluated the hole energy for different values of A and J. The re- 
sults (in MeV) are reported in table [T] For sake of simplicity we have assumed 

M A = A(m + n) (104) 
and the chemical potential /x is chosen as usual as /i = —8 MeV. 
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A 


J = 1/2 


J = 3/2 


J = 5/2 


J = 7/2 


J = 9/2 


12 


-9.5 


-8 








24 


-10.5 


-9 


-8 






40 


-12.5 


-10.5 


-9 


-8 




60 


-15 


-13 


-11.5 


-10 


-8 



Table 1: Energy levels of the relativistic "shell model" for different A and J 



The above example shows how our formalism works. To our knowledge 
the approach presented in this paper is beyond the usual relativistic shell 
calculations, since the usual ways to afford relativity (see our ref. |3j and 
the many references quoted therein) mainly concern QHD (Quantum-Hadro- 
Dynamics) inspired models with a space-dependent mass term that explicitly 
breaks Poincare invariance. This flaw is obviously not obnoxious when heavy 
nuclei are concerned, it has no future however when handling a nucleon as a 
3-quark system. 

7 Conclusion and outlook 

In the present paper we have shown how a relativistic theory of the nucleus 
can be constructed still preserving the main features of the shell model. 
In our approach in fact a shell model like equation has been constructed, 
admittedly in a well defined reference frame, but we have also built up all 
the formalism needed to boost the results to any other frame of reference, thus 
reconstructing Lorentz and Poincare invariance: this is a by far non-trivial 
achievement, since in the traditional nuclear physics translational invariance 
is broken from the very beginning by the shell model even in a non-relativistic 
scheme. Of course the above is particularly suitable for small systems, since 
recoil and center-of-mass motion is fully accounted for. This goes clearly 
beyond the approaches based on translationally invariant systems [21 03 • 

Of course some approximations can be needed in practical calculations, 
and mainly we introduce a "shell-model potential" which is thought to ap- 
proximate the mass operator. Again, exactly as described in sec. |2]we can 
use the same potential to describe particle and hole dynamics, but still the 
same disease survives, since in principle the mass operator in the "hole" and 
"particle" kinematics are intrinsically different. Thus we can use the same 
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potential as a starting point, but then different perturbative expansions are 
required, as shown in sec. 03 

The key issue of the paper is the definition of the "shell model wave 
function": we systematically use quotation marks in referring this quantity 
because it is not at all a wave function, but is defined, instead, as the ex- 
pectation value between physical states (containing any kind of particles, 
namely nucleons, mesons and antinucleons) of some field operators. Thus 
these quantities, referred to in the above as ip and Xi maintain the formal 
analogy with the the true nuclear wave functions of the non-relativistic shell 
model, but contain a much more involved dynamics, since as many mesons 
and antinucleons as possible are allowed to appear, and hence accounted for, 
inside ip and Xi t ne only constraint being a variation of the baryonic number 
of =pl. 

Thus we can apply the formalism developed so far to any relativistic 
system, not necessarily to nuclei, but also (as obvious) to nucleons, where 
the "shell model wave function" (still with quotation marks!) can be regarded 
as the analogous of the quark wave functions in the constituent quark model, 
not disregarding, however, the parton content of the constituent quark, which 
is a composite object built up on current quarks, antiquaks and gluons. 

Our formalism can also be regarded as a theoretical ground for the con- 
stituent quark model and at the same time shows its limitations: in fact, as 
shown above, we can derive from the ip and x the static properties of the 
nucleus or of the nucleon, but the response functions (in the nuclear case) 
or the 7 or lepton interaction with a nucleon require a more detailed study, 
since the degrees of freedom embodied in the "shell model wave function" 
require to be explicitly dealt with. We plan in a successive work to explore 
the dynamical properties of a relativistic complex but finite system. 
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